Hall effects in Bose-Einstein condensates in a rotating optical lattice 
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Using the Kubo formalism, we demonstrate fractional quantum Hall features in a rotating Bose- 
Einstein condensate in a co-rotating two-dimensional optical lattice. The co-rotating lattice and trap 
potential allow for an effective magnetic field and compensation of the centrifugal potential. Frac- 
tional quantum Hall features are seen for the single-particle system and for few strongly interacting 
many-particle systems. 



I. INTRODUCTION 

Rotating Bose-Einstein condensates (BECs) have been 
a subject of great theoretical and experimental inter- 
est over the last few years. Starting from the quan- 
tum engineering of a single vortex rotating conden- 
sates have been used to understand exotic phenomena 
such as the formation of Abrikosov vortex lattices 0, 0| 
and the BCS-BEC crossover Parallels have been 

drawn between a rotating condensate and electrons in 
the presence of a magnetic field and, in particular, the 
fractional quantum Hall effect (FQHE) has been pre- 
dicted for a two-dimensional (2D) condensate rotating 
at a frequency matching that of the confining harmonic 
trap [1, 0, 0, B B • However, the strongly correlated 
FQHE regime has eluded experimentalists in cold quan- 
tum gases due to two problems: it is difficult to con- 
fine condensates at rotation speeds matching the trap- 
ping frequency and vortex shears destroy condensates at 
high rotation. A potential solution is the use of a 2D 
lattice. Introducing a co-rotating optical lattice in the 
tight-binding regime, in which particles on a lattice site 
can only tunnel to adjacent sites, provides strong confine- 
ment and enhances interactions to enable entry into the 
strongly-correlated regime. Similar systems have been 
experimentally demonstrated outside the tight binding 
regime In this paper, we present a direct mapping 



between the angular velocity of a rotating condensate f2 
and the magnetic field as characterized by a in standard 
condensed matter literature on the quantum Hall effect 
(e.g. [l^). Connections have previously been made to 
the FQHE for cold atoms in a lattice in the presence of 
an effective magnetic field [H, or induced tunneling 
loop s p^ . Further, in a recent paper, Umucalilar et. al. 
[l6| presented the phase diagram for bosons in an optical 
lattice in the presence of an effective magnetic field. 

We study a BEC in a rotating 2D optical lattice using 
a Bose-Hubbard Hamiltonian modified by the rotation. 
Observables are computed using exact diagonalization 
with box boundary conditions. There are two advantages 
to using box boundary conditions over periodic boundary 
conditions. First, periodic boundary conditions are suit- 
able only for values of a in a narrow region around ratio- 
nal values while box boundary conditions can be used 
to study the system response for both rational and irra- 
tional values of a. Second, non-periodic elements (such as 
a trapping potential or a lattice tilt) can be introduced 



easily. In the case of box boundary conditions, these 
advantages come with the twin costs of non-negligible 
boundary effects and limited access as the study of many- 
particle systems via the exact-diagonalization method 
quickly becomes intractable with either an increase in 
lattice size or an increase in the number of particles. 

We use the Kubo formalism to describe the system's 
current and density responses to a perturbative potential 
gradient. Using a high-frequency perturbation to over- 
come finite-size effects, we observe FQHE features in a 
single particle system. In particular, the system demon- 
strates plateaus in the transverse resistivity concurrent 
with dips in the diagonal resistivity for fractional values 
of a. At these same values, in concordance with Jaksch 
et. al. [ij, the site number density is modulated with 
a periodicity 1/a. Numerical results are also presented 
for small many-particle systems. In direct analogy with 
the classical Hall effect, a pileup of particles due to the 
Coriolis force is seen along the transverse direction. 

The need for theoretical methods to study the strongly- 
correlated FQHE regime for bosons in a rotating optical 
lattice is urgent, as experimental capabilities to realize 
such systems are rapidly coming to bear. Two main 
characteristics of this regime are site number densities 
of order unity due to strong interactions and filling fac- 
tors (particles per vortex) of order unity needed for the 
creation of composite particles necessary to observe quan- 
tum Hall phenomena. Accordingly, the first experimental 
requirement is for an optical lattice in the tight bind- 
ing regime, traditionally with lattice spacing d ~ 0.5/im. 
The second requirement is for the energy associated with 
rotation to be of the order of the lattice recoil energy or, 
equivalently, for the associated Larmor radius to be of 
the same order as the lattice spacin g. B oth requirements 
have been separately satisfied pTi. Il8l|. An immediate 
benefit of such experiments is the measurement of the 
equivalent of the flux quantum (e^/?i ) constant for mass 
transport. 

This paper is structured as follows: Section II presents 
a derivation of the modified Bose-Hubbard Hamiltonian 
used to study the lattice system along with a discussion of 
current operators. This Hamiltonian is closely connected 
to that traditionally used to study Bloch electrons in a 
magnetic field, but a detailed derivation is useful given 
the new context of cold gases. Section HI describes the 
Kubo formalism used to study the linear response of the 
system. Section IV presents single-particle results for 
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large lattices. Section V contains response characteris- 
tics for small many-particle systems. The last section 
discusses experimental implications and the future out- 
look for this problem. 



II. BOSE-HUBBARD HAMILTONIAN IN 
ROTATING FRAME COORDINATES 

The derivation of a modified Bose-Hubbard Hamilto- 
nian using the symmetric gauge is presented in this sec- 
tion. The angular velocity ft is mapped onto the param- 
eter a used commonly in quantum Hall literature. The 
comparison of this system with that of Bloch electrons 
in a magnetic field is explored by looking at the single- 
particle energy spectrum. 



A. Derivation 

The system to be described is a cloud of a fixed number 
of bosons rotating with an angular velocity about the 
z-axis. This cloud is trapped in a 2D optical lattice co- 
rotating with the same angular velocity in the presence of 
an additional, superimposed two-dimensional harmonic 
trapping potential of frequency to. For a non-rotating 
system {D, = 0), the Hamiltonian Hq has components 
corresponding to the kinetic energy, the lattice and har- 
monic trapping potentials, and the energy due to interac- 
tion between particles. The effect of rotation is included 
by using time- independent rotating-frame coordinates by 
means of the transformation, H = Hq — J dx$triLj,$ |20j 
where $ is a bosonic annihilation field operator describ- 
ing the atoms and is the angular momentum operator. 
The Hamiltonian can then be written in rotating frame 
coordinates as. 



H 



I 



2m 



-hF'^*(x) + F*(x) + f $t$ -QLA<i> ,(1) 



where m is the mass of a single particle and g is the 
coupling constant for repulsive two-body scattering via 
a contact interaction. In this paper, we use a square 
lattice potential described by V^^*^{x) = Vb(sin^(7rx/c?) -I- 
sin^ {ny / d)) . Finally, the trapping potential is y*(x) = 
muP'r^ 12 with r = |x|. Equation ([1]) can be rewritten as, 



R = 



(2) 



Here, 11 = —ih\7 + 7nA(x) is the covariant momentum, 
and A(x) = f2 x x is the equivalent of a magnetic vector 
potential stemming from the rotation. 



The field operator $ can be expanded in several ways. 
One common expansion for the stationary lattice prob- 
lem uses Wannier orbitals Wg{x — Xj), where the sites 
are indexed by i and the bands by I [22l |. If the energy 
separation between the lowest Bloch band and the first 
excited band is large compared to the interaction energy 
and the angular velocity is low {Tifl ^ O.QlEji), then a 
good description is obtained by retaining only Wannier 
orbitals constructed from the lowest Bloch band, i. e., 
/ = 0. With this approximation, the phase description 
of the single-particle wavcfunction is flat within a partic- 
ular lattice site with sharp gradients at site boundaries. 
However, for larger angular velocities {h^l ^ Q.lEfj), the 
riL^-term mixes in higher bands to a non-negligible ex- 
tent. The primary effect of this mixing is to modify the 
phase structure within sites. A modified Wannier basis 
given by 



Wr{^ - Xi) = cxp (^-ij I A(x') • dx'^ I^O(x - Xi) 

(3) 

captures some of this effect by making the azimuthal 
phase gradient within a site proportional to Q. The lower 
limit in the integral is chosen to coincide with the site 
center. This choice ensures that at x = Xi the Wannier 
orbital is real, i. e., Wi?(0) = Ws{0)- A path of integra- 
tion needs to be chosen such that the basis set defined by 
Wii(x — Xi) satisfies orthonormality. In addition, calcu- 
lations for the square-lattice problem are greatly simpli- 
fied if the choice of path allows for each two-dimensional 
Wannier orbital to be decoupled into a product of one- 
dimensional Wannier orbitals. One particular choice of a 
path that meets both criteria is along straight lines paral- 
lel to the lattice axes. Note that this path is not unique 
and decoupling is not required. A quantitative assess- 
ment using imaginary time propagation techniques (23j 
shows that the modified Wannier basis set Wfl(x — Xj) 
describes the phase gradient within a site better than 
the regular Wannier basis Ws(x — Xj) and captures the 
pertinent parts of the Hilbert space for our discussion. 

Using the modified Wannier basis Wi^(x — Xi), the field 
operator can now be expressed as 



l>(x) =^a,W^i?(x-Xi), 



(4) 



where a,; is a site-specific annihilation operator. A modi- 
fied Bose-Hubbard Hamiltonian is obtained by substitut- 
ing Eq. Q into Eq. Q, 



^1 ialaje-''^'^ + a,a]e"^^^) 



+ -^n.i{ni - 1) 



(5) 



where indicates the sum over nearest- neighbor site 
pairs and fii = a|ai is the number operator. The dcfi- 
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FIG. 1: Overlap integrals Ai and A2 as a function of 
lattice depth for a standing-wave optical lattice described 
by l/('''''(x) = Vo(sin2(7rx/d) + sin'^ {ny / d)) . The lat- 
tice depth is given in units of the recoil energy Er — 
Tt?-k^ /2md? . A-i and A2 are in units of d? . In the tight- 
binding regime {Vo/Er > 5), Ai ~ -0.273 exp (-Fo" ®^*') and 
A2 ~ 0.368 exp {-V^-^^''). 



nitions for </)y, e, Ai and A2 follow. The phase for the 
hopping term is 

m f^' , , mO 

(Pi] = Y ) ■ ax = —{x,yj - Xjyi) . (6) 

Here, the angular velocity CI is expressed in units of En/h 
where En is the recoil energy associated with lattice spac- 
ing d. The parameters t and e are identical to the hop- 
ping and onsite zero-point energies associated with the 
standard Bose-Hubbard model [24| and are obtained by 
evaluating the integrals 



t 



J d^W^{K - xi) (-|^V2 + y('»*)(x)) Ws{k - Xj) , 



(7) 



J dxVFl(x-Xi) (-^V2 + y(''^*)(x)^ M^s(x-Xi). 



(8) 



The modifications to these terms due to rotation are pro- 
portional to (f2^— w^) and to two new overlap parameters 
given by 

Ai = [ dxW^{x-x,) ix-x, f Wsix-xj), (9) 



dxW^ix ~ X,) (x - x^f Wsix - X,) ,(10) 



where Ws{x — Xi) is a one-dimensional Wannier func- 
tion. There is an additional factor of two in A2 because 
of identical onsite overlaps along the x- and y- directions. 



Changes in the lattice potential affect both new parame- 
ters (Fig.[T]). The last term in the Hamiltonian describes 
the onsite-interaction energy, and for an s-wave scatter- 
ing length as [li,[2|. 



U = 



J dx|iys(x-xi)|' 



(11) 



The interaction term as described in Eq. (jll[) is included 
for completeness and is used only implicitly in this paper. 
The single-particle discussions are trivially independent 
of U, while the many-particle hardcore-boson analyses 
implement the U ^ 00 condition using the two-state 
approximation. 

There are two other approximations implicit in our 
approach. The first is the tight-binding approximation 
whereby only hopping between adjacent sites is consid- 
ered. This approximation becomes valid when Vq ex- 
ceeds 5E]i. Our calculations are well in the tight-binding 
regime with Vo = lO-E^. The second approximation is 
the use of infinite-lattice Wannier functions for a finite 
lattice. Due to this, edge effects are not accounted for 
correctly though the approximation gets better with in- 
crease in lattice size. 

We implement this approach by constructing the 
Hamiltonian using a truncated Fock-number basis for the 
desired number of sites and diagonalizing it numerically. 
Note that we implicitly introduce infinite potential walls 
around the lattice by spatially restricting particles to a 
limited number of sites. As will be discussed in Sec. IIIH 
we minimize the effect of the box boundary conditions in 
the linear response by modulating the perturbation with 
a period much smaller than the time scales associated 
with particles tunneling from one lattice site to the next. 

Note that the modified Bose-Hubbard Hamiltonian in 
Eq. ([5]) is obtained using the symmetric gauge for the 
vector potential A(x). An expression equivalent to that 
used by Jaksch et. al. [l^ can be obtained using the Lan- 
dau gauge with Al(x) = flyi and setting the trapping 
frequency equal to the angular velocity (lu = f2). The 
two different Hamiltonians can be connected using the 
transformation Hl = PlHP£^ where 



Pl 



exp 



, mf2 



x'dy' 



(12) 



B. Mapping angular velocity onto a 

Several two-dimensional problems are characterized by 
multiply-connected domains where singularities in the 
topology are typically due to quantized magnetic flux 
lines (e. g. Aharonov-Bohm effect [l^l) or strongly re- 
pulsive particles (e. g. quantum Hall effect [13, [23). In 
this context, it is useful to introduce a winding number 
2T:a that describes the phase picked up by a particle when 
it goes around such a singularity. Inaccessible regions in 
the topology can also be created by means of a suitable 
potential. 
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FIG. 2: (Color online) Schematic for a particle going around a 
plaquette in a square lattice. Regions shaded dark correspond 
to lattice sites and the light regions indicate peaks in the 
lattice potential. Using Eq. ((6}, it can be shown that the 
particle picks up a phase of 2-Ka — 2md?Q/h as it returns to 
its original position as marked by the solid circle. If the path 
of the particle encloses P plaquettes then the phase picked up 
is 2naP. 



Consider the lattice potential shown in Fig. [2] The 
light shaded regions correspond to peaks in the lattice 
potential that are inaccessible to particles in the tight 
binding regime. The phase accumulated by a particle 
adiabatically going around one such simply connected in- 
accessible region (a plaquette) can be calculated by first 
breaking the loop into four parts as indicated. For each 
part, the phase change associated with destroying a par- 
ticle at a site and creating it in a neighboring site is given 
by Eq. ([6]) . This phase is identical to that associated with 
the hopping term in the Hamiltonian [Eq. ([5|)]. The rela- 
tionship between the angular velocity and a is obtained 
by summing the contributions and is given by 

Henceforth, we will use a to characterize the rotation 
frequency instead of the angular velocity fl in order to 
maintain the connection with quantum Hall literature. 

C. Single-particle energy spectra 
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FIG. 3: Single particle energy spectra as a function of the 
rescaled angular velocity a for a 40 x 40 lattice. Darker shad- 
ing indicates greater density of states. The origin of the y-axis 
has been shifted to coincide with the onsite energy. 



context has two consequences for the energy spectrum: 
first, the spectrum becomes periodic as a function of a 
with a periodicity Aa = 1 and second, the spectrum is 
symmetric about a = 0.5 and takes on the shape of the 
Hofstadter butterfly [2^ — originally used to describe the 
energy spectra for an electron in a periodic potential in 
the presence of a magnetic field. For oj ^ fl, both the 
symmetry of the energy spectra about a = 0.5 and the 
periodicity are disrupted as the entire spectrum shifts up 
or down as a function of (il^ — w^). As shown by Ana- 
lytis et. al. [s^l , the fractal nature of the spectra becomes 
increasingly well-defined as the size of the lattice under 
consideration grows. 

The grayscale in Fig. [3] describes the density of states 
and the finite nature of the lattice manifests itself in the 
sparse energy levels between bands. For a ^ 1, the 
lowest bands are linearly proportional to a leading one 
to draw comparisons to the Landau energy spectra for a 
single particle in a 2D harmonic oscillator. The Landau 
energy levels are given by En = {n -\- \/2)fujJc where n is 
an integer and u>c is the cyclotron frequency. For small a 
in the spectra shown in Fig. [3l the slopes are not similarly 
proportional to the band index, e. g., for the five lowest 
bands, the slopes are ~ 6, 17, 27, 35 and 43. 



Having derived the Hamiltonian, a useful crosscheck 
is the comparison of the single-particle energy spectrum 
with that for a Bloch electron in the presence of a mag- 
netic field. The energy spectrum for a single particle in 
a 40 X 40 lattice is plotted as a function of a in Fig. [3l 
The energy contribution due to the centrifugal force is 
eliminated by setting the trapping frequency equal to the 
angular velocity (a; = r2). This is identical to the con- 
dition necessary to reach the highly degenerate lowest 
Landau level (LLL) for the same problem in the absence 
of a lattice. Applying the LLL condition in the lattice 



III. KUBO LINEAR RESPONSE 

The Hall effect describes the longitudinal and trans- 
verse transport responses of a two-dimensional electron 
gas in the presence of a magnetic field to an applied elec- 
tric potential gradient. The mapping between magnetic 
flux density for the 2D electron problem and angular ve- 
locity for a rotating gas is valid when the latter problem 
is formulated in rotating frame coordinates. Accordingly, 
a potential gradient is introduced in the rotating frame 
by linearly modifying the lattice onsite energy along the 



5 



direction of the tilt. This section is divided into two 
parts. The first part lays out two quantities useful for 
studying particle transport in this system: the inter-site 
current operator and onsite density. The second part 
briefly sketches a derivation of the Kubo formula used to 
study the linear response of this system. 



A. Current and density operators 

The single particle current in the rotating frame, equiv- 
alent to the Nother current associated with local phase 
changes of the wavefunction, is obtained by using the 
mass continuity equation for an infinitesimal volume. 
The current operator is then realized by quantizing the 
field using Eq. (jl]) and is given by, 



i^(x) 



1 
2m 



$t(x) ( -V - mA(x) ) $(x 



V - mA(x) $(x) $(x) 



(14) 



The connection to the current operator for the station- 
ary lattice is made by using the transformation $(x) = 
Pl>s(x)P"^ where. 



P = exp 



• dx' 



Substituting for <l>(x) in Eq. HH) yields, 



J^(x) = — /m ($l.(x)V$s(x) 



(15) 



(16) 



P(j^(x)+A(x)$t(x)$(x)) P-I, (17) 



where the second term in Eq. (|17[) is needed for conser- 
vation of current [3ll |. We obtain the current from site 
i into site j by integrating the Nother current [Eq. (fT4|l ] 
across the boundary between the two sites, 

im 



OjCiie 



(18) 



where 



d/2 



~d/2 



dyW^{y)Ws{y) 



Ws[x~-] d.,Ws I X 



(19) 



<2 



For a lattice depth of Vb = lO-Bfl, 7 w 0.094d^ 

The onsite density operator for a unit cell i [Fig. [5] is. 



(20) 



where Bi and B2 are dimensionless overlaps within a unit 
cell for Wannier functions centered on the same site and 
adjacent sites respectively. For our calculations well in 
the tight-binding regime, Vq = IQEj^, Bi w 0.9969 and 
Bo « 1.6591 X 10"^ 




FIG. 4: (Color online) Scheme for creating a perturbative lin- 
ear gradient potential in the system. The lattice potential is 
tilted by an angle Q along the a;-direction. The 2D trapping 
potential is cancelled out by the centrifugal force at = f2 
and has not been shown. In order to reduce the effect of the 
implicit infinite boundary potential walls in our calculations, 
we consider an AC perturbation where the tilt angle is mod- 
ulated by a frequency — dO/dr. 



B. Kubo formula 

Consider the lattice system sketched in Fig. ID The op- 
tical lattice is co-rotating with the condensate about the 
2;-axis with angular velocity f2. The harmonic trap (not 
shown) frequency is adjusted to a; = f2 such that the cen- 
trifugal force is cancelled. The perturbation is introduced 
by tilting the lattice along the ai-axis in the rotating frame 
and is modulated by a frequency ly to induce sloshing. An 
AC perturbation is switched on at time r = 0. A com- 
mon mathematical trick to simultaneously extract both 
quadrature components of the linear response of the sys- 
tem is to use a complex perturbation. In this case, the 
sine and cosine (phase-shifted) components of V will go 
through and recombine to give a exp(wr) factor in the 
final result. Accordingly, the perturbation is written as. 



V{t) = AeiT)e"'^X = Ae(T)e"^^ ^ a 



(21) 



where A is the strength of the perturbation, Xi is the a;- 
coordinate of site i, and Q{t) is the Heaviside function. 
The effect of the implicit infinite boundaries is mitigated 
by making the time scales associated with the sloshing 
small compared to that associated with hopping from 
one site to the next, i. e., hv >> t. A brief sketch of 
the derivation for the change in the expectation value 
of an observable Y due to the perturbation follows (see 
Ref. [H, [111 for detailed discussions). 

The density matrix in the interaction picture {t) can 
be broken into a time-independent part and the change 
Ap^(r) due to the perturbation. 



(22) 
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where the superscript / marks quantities in the interac- 
tion picture. The time independent part po corresponds 
to the density matrix for the unperturbed system. Re- 
taining the first order terms in the Liouville equation of 
motion for the density matrix provides an expression for 
the second term in Eq. 



~i;(t-t') 



dr'. 



(23) 



Here 77 is used to fix the boundary conditions and we 
take the hmit 77 — > 0+ at the end of the calculation. The 
expectation value of Y is, 

(y(r)) = Tr {y{t)p{t)} = Tr {y\t)p'{t)} . (24) 

The expectation value of the response to the perturbation 
is, 



(AYir)) =Tr[r{T)Ap'{r)} . 



(25) 



At low temperatures, the only contribution to the trace 
comes from the ground state, i.e., po « |'0o)('0o|- The 
final expression for the expectation value of the response 
is obtained using this approximation and substituting for 

Ap^(T) [Eq. mi 



{AY{r)) 



■E 

ri>0 



(V^0|^|^n)(^n|^|V'0) 



{un - ujo + v) - iri 



(26) 



where 



{^0 



X 



2pn) are energy eigenstates. Note that n > 0, i.e., 
■00) = , because the unperturbed ground state 



is symmetric about the y-axis while X is not. For the 
purposes of this paper, the hnear response [Eq. ((26|) ] is 
evaluated at limits t — > 00, 77 — > 0-1- and the prefactor 
exp(ii/T) is excluded from results shown. 



IV. SINGLE-PARTICLE RESPONSE 

This section presents numerical results for the linear 
transport response of a single particle in a rotating 40 x 
40 lattice. The system is subjected to a perturbation 
modulated at frequency v ~ En/fi. The linear response 
is characterized in terms of the change in end currents 
and the sample averaged resistivity. 

The ideal way to study the current and voltage charac- 
teristics of the system would be to connect it to reservoirs 
and compute currents between the system and reservoirs, 
as is done in studying open quantum systems. However, 
this approach becomes numerically intractable for sys- 
tems of size > 4 — 6 sites [13] ■ For an isolated system. 




Field gradient 



FIG. 5: (Color online) Top view schematic depicting end cur- 
rents for a 8 X 8 lattice. The arrows crossing the solid lines 
mark the longitudinal end currents {J^) while the arrows 
crossing the dashed lines indicate the transverse end currents 
{Jy)- The lattice is tilted to the right. 



proxies for the in and out current response of the system 
are the end currents — the current response of the system 
very close to the boundaries of the system. Note that this 
is true only for a linear response study. The operators for 
the end currents are obtained by summing current oper- 
ators across end links as shown in Fig. [5] End currents 
along each direction are added on either side of the lat- 
tice in order to capture only additional currents due to 
the perturbation. The underlying currents/circulation of 
the system due to the rotation have been described else- 
where [35j . 

The expectation values for the end currents along the 
longitudinal and transverse directions are plotted as a 
function of the winding rate a [Eq. in Fig. [51 The 

longitudinal end current displays well-defined peaks at 
fractional values of a. At these fractional values of 
a = p/q (where p,q are integers), the energy spectra 
[Fig. in] breaks up into exactly q bands [2§|. For low 
perturbation frequencies {hv ~ t), the denominator in 
Eq. p6|) is very small for nearly degenerate states within 
the same band and the linear response is, in general, 
large. The system described in this paper has implicit 
infinite potential walls and the perturbation frequency 
is far off resonance {hv = Eji ^ 50t) in order to elim- 
inate Bloch oscillations. Therefore, the denominator in 
Eq. 126]) does not become resonant for any value of a. 
The peaks appear due to bigger off-diagonal current ma- 
trix elements at fractional values of a and are small due 
to their non-resonant character. For high frequencies, the 
height of the peaks goes as 1/v. 

The plot of the transverse current [Fig. [6l^b)] displays 
peaks/dips at the same values of a. The transverse end 
current is antisymmetric about a = 0.5. To understand 
this, consider a value of a = 1 — /3. The corresponding 
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As the size of the lattice under consideration gets big- 
ger, the peak structure in Fig. [6] becomes more well- 
defined in two ways. First, the peaks become narrower 
as they get centered closer to exact fractional values of a, 
and second, more peaks appear at other fractional values 
of OL. Both of these effects correspond to better resolution 
of the fractal nature of the energy spectra with bigger 
lattice size. The height of the peaks, however, decreases 
exponentially with lattice size. For example, consider the 
current plotted in Fig. [H The height of the central peak 
goes as ^ 1.7exp(— 0.6L) where L is the number of sites 
along a side of the lattice. 

A spatial average of the current response across the 
system smoothens out the peaks in Fig. [6l The conduc- 
tivity tensor describes the response of the sample aver- 
aged current. If the perturbation is along the x-direction, 
the conductivity tensor elements are 



FIG. 6: Expectation value of end currents, (a) ^3^) (along 
the x-direction) and (b) (Jjf) (along the y direction), as a 
function of q for a single particle in a 40 x 40 lattice subject 
to a linear ramp perturbation of amplitude A modulated with 
a frequency v = Ea/h. 
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FIG. 7: (a) Diagonal and (b) transverse resistivity [Eq. (1281) ] 
as a function of angular velocity for a single particle in a 40 x 
40 lattice subject to a linear-ramp perturbation of amplitude 
E modulated with a frequency v = En/h. The shape of 
the diagonal resistivity plot resembles the cross-section of the 
Mott Insulator lobe at /i/(7 = 0.5 as seen in Ref. [l^ . 



angular velocity is 17 = {nh/M(P){l - /3) [Eq. dTS])]. The 
phase picked up by a particle going around a plaquette 
is 27ra = {2tt - 27r/3) ee -2tt(3 [Fig. [2]. The latter phase 
winding can equivalently be created by rotation in the op- 
posite direction with angular velocity = ~{'Kh/Md?)(3 
for which the Coriolis force v y. Vl) is in the opposite 
direction. 



A 



(27) 



where AJ^ indicates linear response of the sample av- 
eraged current along the /i direction (i. e., the total re- 
sponse of all current operators for links along the ^ di- 
rection) . The resistivity tensor elements are derived from 
the conductivity tensor using. 



(28) 



The sample-averaged longitudinal and transverse resis- 
tivities are shown as a function of a in Fig. [71 The plot 
of the longitudinal resistivity has dips at all fractional 
values of a, though the dips are now seen only around 
prominent fractions such as a = 1/2, 1/3, 2/3, . . .. These 
fractions arc the most common in the sense that for a 
given range of integers, these fractions can be constructed 
in the most number of ways. The plot of the transverse 
resistivity shows plateaus at values of a corresponding 
to these dips. Both features are signatures of the FQHE 
seen in a 2D electron gas. This appearance of a many- 
particle effect in a single-particle system is intriguing. A 
tentative explanation is given by considering the effect 
of the optical lattice. In a 2D electron gas, the com- 
bined effect of the magnetic field and the Coulombic in- 
teraction is to arrange the electrons into a lattice. For a 
filling factor of one, the electrons fill the lowest Landau 
level forming a hexagonal lattice in the nearest neighbor 
approximation. The lattice spacing is 2^^tt /Hb^ where 
Ib = \/Ti/ eBi_ is the magnetic length determined by a 
magnetic field B±^ [l^. In addition, the two-dimensional 
geometry in which electrons cannot cross each other leads 
to a change in phase equal to 27ra each time one electron 
circles another. These effects are reproduced when a lat- 
tice is introduced in such a way that the particle picks 
up a phase of 27ra going around a plaquette [Fig. [2]. In 
Fig. [51 square periodic density structures are seen for 
the single-particle case at certain values of a in where 
the expectation value of the ground-state site number 
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FIG. 8: Unperturbed ground state site number-density distributions for a single particle in a 40 x 40 lattice for (a) a — 1/4, 
(b) Q = 1/3, (c) a = 2/5, (d) a = 1/2, (e) a — 3/5, (f) a = 2/3, (g) a = it/11, and (h) a — 1. For simple fractions such 
as a = 1/4 (a) or Q = 1/3 (b), the site number-density distribution has peaks separated by 4 and 3 sites respectively. For a 
fractional value such as a = 2/5 (c), the number distribution has periodically arranged rings with centers separated by 5 sites 
[Fig. HJc)]. The density distributions for any value of a are the same as those for 1 — a (compare (b) and (f) or (c) and (e)). 
For non rational values, e. g. a = tt/II (g), the periodicity is complicated, if not destroyed. The 2D Gaussian-like envelope seen 
in all the subplots are due to 2D box infinite wall conditions. 



density for the unperturbed system has been defined as 
(fii) = {a\ai). At values of a = p/q {{p,q} G integers) 
corresponding to dips in the longitudinal resistivity, {fn) 
has a periodicity q — nh/m(Pil. This stems from the 
periodicity due to the hopping term in the Hamiltonian 
[Eq. ([2])] (see also Refs. [ij, O [T^). Bragg scattering 



is a promising probe for such structures [13|. Note that 
the periodicity goes as l/D, and not as as might 

be expected by direct comparison with the electron gas 
system. 



As discussed in Ref. [ij], for a ^ 1 in an infinite sys- 
tem, the length scale of the wavefunctions is much larger 
than the lattice spacing and in this continuum limit, the 
ground state of the system is the one-half Laughlin state. 
In addition, the site number density distributions for a 
and 1 — a are identical (e. g., Figs.[8l^b) and (f)). There- 
fore, the number density distribution for a = 1 is the 
same as that for a = and corresponds to a system 
without rotation or lattice. The concentration of parti- 
cles at the center [Fig. [8l[h)] is due to the 2D infinite box 
potential. 



V. MANY-PARTICLE RESPONSE 

The introduction of more hard-core bosons to the two- 
dimensional system described earlier adds an additional 
degree of freedom to the problem. In the single-particle 
system, circling a plaquette added a phase of 2T:a. This 
is still true in the many-particle system but there is an 



additional phase when two particles are exchanged [36l |. 
i. e., 4'(xi,X2) = exp(i27ra')\E'(x2, Xi), where a' is a di- 
mensionless winding rate similar to a. This substantially 
complicates the problem. This section extends the earlier 
analysis to many particle systems using numerical results 
for small systems. 

The energy spectrum for two hardcore bosons in a 8 x 8 
lattice is plotted in Fig. [9l The overall butterfly outline 
seen in Fig. [3] is preserved. For N particles, the total 
energy band width defined by the maximum energy dif- 
ference at a = 0, is AE^ax = 8tN. The gray shading 
describes the density of states which is marked by de- 
generacy at energy E = 0. It is difficult to delineate 
a band structure due to finite lattice size. However, at 
the most distinct regions {a ^ 0.5), there appear to be 
three bands as opposed to two seen in Fig. [3l Diagonal- 
ization of larger many-particle systems quickly becomes 
intractable due to the exponential scaling of the Hilbert 
space dimension with particle number. 

A useful tool for understanding the structure of a 
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FIG. 9: Energy spectrum for two particles in a 8 x 8 lattice. 
The gray shading marks the density of states on a logarithmic 
scale. The trapping frequency uj is set equal to the angular 
velocity fi. 



many-body state is the long range order, which, for a 
lattice system, can be defined by 



(29) 



The hardcore natm'c ([/—!• oo) of the particles manifests 
itself as an anti-correlation envelope in the ground state 
seen most clearly for a = [Fig. [TO^a)]. The periodicity 
in the Hamiltonian for rational a discussed earlier is seen 
in the long range order of the ground state plotted for 
a = 1/2 and a = 1/3 [Figs.[10i;c) and[10te)]. The Fourier 
transforms of the long range order [Figs.fTOFd) andfTOTf)] 
show peaks at g = l/2c? and q = l/Sd corresponding to 
periodicities of two and three lattice sites respectively. 
Interestingly, every second site is correlated for a ~ 1/2 
while every third site is anti-correlated for a = 1/3. 

The features of the end current response [Fig. [S] 
as a function of a for two particles in a 8 x 8 lat- 
tice are altered considerably due to finite-size effects 
[Fig.[n]- The longitudinal {{J^)) and transverse {{Jy)) 
end currents display similar features close to a = 
1/2, 1/3, 2/3, 1/4, 3/4 . . . though the transverse end cur- 
rents are antisymmetric about a ~ 0.5. Both single- 
particle and two-particle end current responses are sim- 
ilar with a few additional peaks in the latter. In the 
single-particle analysis, the first distinct peaks {a ~ 
1/2, 1/3, 2/3) start emerging for lattice sizes > 10 x 10, a 
size just beyond our numerical methods for two particles. 

The resistivity [Eq. p8|) ] scaled by the number of par- 
ticles as a function of a is plotted in Fig. [121 In the low 
filling limit (particles/site < 0.1), the effect of the interac- 
tion between particles is to enhance the scaled resistivity 
or, equivalently, lower the sample-averaged conductivity 
per particle. This is consistent with earlier findings [35| . 
where increasing interaction reduced the current per par- 
ticle. The two particle resistivity also shows weak dips 
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FIG. 10: Ground state long range order [Eq. (|29p ] in position 
(x) and quasi-momentum space (q) for two particles in a 16 x 
16 lattice. The ground states are calculated for a = [(a) and 
(b)], a = 1/2 [(c) and (d)], and a = 1/3 [(e) and (f)]. One of 
the two points Xi is fixed close to the center of rotation while 
Xj is moved from end to the other parallel to one of the edges 
of the lattice. The asymmetry in (a),(c) and (e) is due to Xi 
being closer to one of the infinite walls. 



(inflections) in the longitudinal (transverse) resistivity at 
fractional values of a. Note that this calculation is in the 
very dilute limit, where the physics is largely dominated 
by single-particle effects. 

One feature of the Hall effect is the breaking of number- 
density symmetry along the y-axis despite the perturba- 
tion being along the x-axis. This leads to charge buildup 
and eventually creates a stopping Hall potential. The 
number density asymmetry can be quantified by the dif- 
ference in the number of particles on each half, where the 
lattice is divided in two along the direction of the per- 
turbation. This quantity Ap is plotted as a function of a 
for four particles in a 4 x 4 lattice in Fig. [131 The effect 
of interaction on Ap is examined by comparing the weak 
(U = 0) and strong (U — oo) interaction limits. In the 
strongly interacting regime, the redistribution (as quan- 
tified by Ap) due to the Coriolis force is three orders of 
magnitude greater than that in the non-interacting limit. 
The change in the direction of particle pileup at a = 0.5 
marks the change in the direction of the Coriolis force. 
Since the particles are charge neutral in this system, the 
retarding potential along the transverse direction is cre- 
ated by strong repulsive interaction between particles. 
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FIG. 11: End current response per particle for one (dashed, 
n = 1) and two (solid, n = 2) particles in a 8 x 8 lattice along 
the (a) longitudinal direction — ' '^^^ 
transverse direction — {Jy)/n. 
lated at a frequency v = Er/Ti. 



{J^ )/n, and (b) along the 
The perturbation is modu- 




FIG. 12: Transverse and diagonal resistivity scaled by the 
number of particles n as a function of a for one (dashed) 
and two (solid) particles in a 8 x 8 lattice. In the low filling 
(particles per lattice site) limit, the effect of the interaction is 
to decrease the conductivity per particle. The perturbation 
is modulated at a frequency u — Eu/h. 



VI. CONCLUSION 
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FIG. 13: Differential density response for four particles in 
a 4 X 4 lattice in the strongly repulsive (solid) and non- 
interacting (dashed) limits. The two dashed lines describe 
the same limit but are displayed at different magnifications. 
In the strongly interacting limit, the redistribution due to 
the Coriolis force for this particular system is three orders of 
magnitude greater than that in the non-interacting limit. 



tional quantum Hall features. Density redistribution in 
small, strongly-correlated many-particle systems shows 
the equivalent of the classical effect and the mapping be- 
tween the Coriolis and Lorentz forces. However, larger 
many-particle systems need to be considered to find the 
FQHE as described by dips (plateaus) in the diagonal 
(transverse) elements of the conductivity tensor as a func- 
tion of a. Tilting the rotating optical lattice in experi- 
ments such as |ll| opens up this topic in two different 
ways. In the weakly interacting limit, the system can 
be probed to study the Coriolis force in superfluid sys- 
The system described by Tung et. 



tems. The system described by Tung et. al. [ll| obtains 
an ensemble of stacked, identical two-dimensional layers. 
If the layers are stacked along the z-axis and the lattice 
is tilted along the a;-axis, the effect of the Coriolis force 
would be to break the two-fold symmetry in the density 
distribution along the j/-axis. An observable to charac- 
terize this is obtained by imaging the density distribution 
along the z-axis using a CCD camera and evaluating Ap, 
the change in density induced by the perturbation. The 
strongly-interacting limit can be achieved by tuning the 
two-body scattering length using a Feshbach resonance. 
In this limit, the ground state can be examined via Bragg 
scattering. Future work will focus on extending results 
obtained here to larger systems by improving the treat- 
ment of boundary conditions. This will be done using a 
combination of periodic boundary conditions and block- 
ing the lattice into 'supercells'. 



This paper considers bosons in a rotating optical lat- 
tice that have a Hamiltonian similar to that for Bloch 
electrons in a magnetic field. The Hall effect in this 
system is probed using linear-response theory and the 
Kubo formula. The single-particle case exhibits frac- 
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